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Abstract 

In the present paper, the temperature distribution in the critical state of hard superconductors is 
investigated in the quasi-stationary approximation. It is shown that the temperature profile can 
be essentially inhomogeneous in the sample, which affects the conditions of initiation of a magnetic 
flux jumps. 
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While dealing with instabilities of the critical state in hard superconductors, the charac- 
ter of the temperature distribution T(x,t) and that of the electromagnetic field E(x,t) are 
of substantial practical interest [1]. This derives from the fact that thermal and magnetic 
destractions of the critical state caused by Joule self-heating are defined by the initial tem- 
perature and electromagnetic field distributions. Hence, the form of the temperature profile 
may noticeably influence the criteria of critical-state stability with respect to jumps in the 
magnetic flux in a superconductor. Earlier (cf., e.g., [2]), in dealing with this problem, it was 
usually assumed that the spatial distribution of temperature and field were either homoge- 
neous or slightly inhomogeneous. However, in reality, physical parameters of superconductors 
may be inhomogeneous along the sample as well as in its cross-sectional plane. Such inho- 
mogeneities can appear due to different physical reasons. First, the vortex structure pinning 
can be inhomogeneous due to the existence of weak bonds in the superconductor. Second, 
inhomogenety of the properties may be caused by their dependence on the magnetic field H. 
Indeed, the field H influences many physical quantities, such as the critical current density 
j c , the differential conductivity a^, and the heat conductivity k. 

In the present paper, the temperature distribution in the critical state of hard superconduc- 
tors is investigated in the quasi-stationary approximation. It is shown that the temperature 
profile can be essentially inhomogeneous, which affects the conditions of initiation of a mag- 
netic flux jumps. 

The evolution of thermal (T) and electromagnetic (E, H) perturbations in superconductors 
is described by a nonlinear heat conduction equation [3], 

^ = V[kVT}+jE, (1) 



a system of Maxwell's equations, 



and a critical-state equation 



* 1 dH 

rotE = (2) 
roW = (3) 



j=l(T,H)+J r (E). (4) 



Here z/ = f(T) is the specific heat, k = k(T) is the thermal conductivity respectively; j c is 
the critical current density and j r is the active current density. 

We use the Bean-London critical state model to describe the j c (T, H) dependence, accord- 
ing to which j c (T) = j [l — a(T — T )] [4], where the parameter a characterizes thermally 
activated weakening of Abrikosov vortex pinning on crystal lattice defects, j is the equilib- 
rium current density, and T is the temperature of the superconductor. 

The j r (E) dependence in the region of sufficiently strong electric fields (E > Ej; where Ej 
is the limit of the linear region of the current- voltage characteristic of the sample [2]) can 

T)C^ <7 n H c 

be approximated by a piecewise-linear function j r m (JfE, where <jf = ~ "^ C2 is the 

Tthc 

effective conductivity in the flux flow regime and n is the viscous coefficient, $o — is the 

2e 

magnetic flux quantum, a n is the conductivity in the normal state, H C2 is the upper critical 
magnetic field. In the region of the weak fields (E < Ef), the function j r (E) is nonlinear. 
This nonlinearity is associated with thermally activated creep of the magnetic flux [5]. 

— * 

Let us consider a superconducting sample placed into an external magnetic field H = 

dH 

(0, 0, H e ) increasing at a constant rate = iJ=const. According to the Maxwell equation 

(2), a vortex electric field E = (0, E e , 0) is present. Here H e is the magnitude of the external 
magnetic field and E e is the magnitude of the back-ground electric field. In accordance with 
the concept of the critical state, the current density and the electric field must be parallel: 
E || j;. The thermal and electromagnetic boundary conditions for the Eqs. (l)-(4) have the 
form 



dT_ 

dx 



+ w [T(0)-T } = 0, T(L)=T , 

x=0 

(5) 
dE 

' = , E(L) = , 

x=0 



dx 



For the plane geometry (Fig.) and the boundary conditions H(0) = H e , H(L) = 0, the 

magnetic field distribution is H(x) = H e (L — x), where L = - — - is the depth of magnetic 

flux penetration into the sample and wo is the coefficient of heat transfer to the cooler at 
the equilibrium temperature T . 
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The condition of applicability of Eqs. (l)-(4) to the description of the dynamics of evolution 
of thermomagnetic perturbations are discussed at length in [2]. 

In the quasi-stationary approximation, terms with time derivatives can be neglected in 
Eqs. (l)-(4). This means that the heat transfer from the sample surface compensates the 
energy dissipation arising in the viscous flow of magnetic flux in the medium with an effective 
conductivity 07. In this approximation, the solution to Eq. (2) has the form 



E = —(L-x). 

c 



(6) 



Upon substituting this expression into Eq. (1) we get an inhomogeneous equation for the 
temperature distribution T(x,t), 

d 2 Q 



dp 2 



-P<d = f(p). 



(7) 



Here we introduced the following dimensionless variables 



f(p) = -[l + rup] 



Jo 

aT n ' 



= 



T-T 
To 



P = 



L — x 



and the dimensionless parameters uo = — — , and, r = ( — = \ , where r characterizes 

cjo \aHL 2 / 
the spatial scale of the temperature profile inhomogeneity in the sample. Solutions to Eq. 

(7) are Airy functions, which can be expressed through Bessel functions of the order 1/3 [6] 
0( P ) = C lP ^K 1/3 (jp 3 / 2 ) + c 2P v% /3 (Jp 3 / 2 ) + 6 (p), (8) 

0o(p) = p^K 1/s (^p 3 / 2 ) [[I + r^pTh/s (Ipf 2 ) d Pl - 



P 1/2 h/, (Ip 3/2 ) [[I + rup^pTK^ (IpT) d Pl , 

where C 1 and C 2 are integration constants, which are determined by the boundary conditions 
to be 



-w LQ(0) + k 



dQ 
dp 



Ci = 0, c 2 = 



p= — 

r 



From the Maxwell equation (2), the temperature inhomogeneity parameter can be expressed 
in the form 

1/3 

(9) 



r 

a= L 



aEl Ht K 
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It is evident that a ~ 1 near the threshold for a flux jump, when 



aHl_ 



1, even under the 



Ht K vL 2 
quasi-stationary heating condition — — << 1; where t K = is the characteristic time of 



K 



the heat conduction problem. 

Let us estimate the maximu 
The solution to Eq. (7) can be represented in the form 



K 



Let us estimate the maximum heating temperature 9 m in the isothermal case w — — > 1. 

Lj 



g(x) = e r 



Po- 



(x Xm) 



(10) 



near the point at which the temperature is a maximum, x = x m (Fig.). 
With solution (10) being approximated near the point x m = — with the help of the 



thermal boundary conditions, the coefficient p can be easily determined to be y— ) m 
and the temperature can be written as 

4 / L x 2 ' 
&(x) = Q m \--lx-- 



Substituting this solution into Eq.(7), the superconductor maximum heating temperature 
due to magnetic flux jumps can be estimated as 



jo + ^—(L-x r , 



H 

ckTq 



(L — x m ) 



1_ 
L 2 



aH 



CK 



(12) 



L — x m ) 



For a typical situation when — << (L — x m ) the estimation for 6 m is 

L CK, 



Jo + — — [L - x m ) 

c 



HL 2 
ckTq 



(L — x m ). 



(13) 



Here, the parameter 7 ~ 1 (for a parabolic temperature profile 7 ~ 8). It is easy to verify 
that for typical values of j = 10 6 A/ cm 2 , H = 10 4 Gs/sek, and L = 0, 01 cm the heating is 
sufficiently low: 6 m << 1. In the case of poor sample cooling, w — 1 — l0erg/(cm 2 sK), the 



G m is 



0m — — ~ 0, 5; 



(14) 



CWqT 

i.e., the heating temperature can be as high as 5T m = T Q m ~ 2K. One can see that 
in the case of poor sample cooling, the heating can be rather noticeable and influence the 
conditions of the thermomagnetic instability of the critical state in the superconductor. 
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FIGURE 

Fig. The distrubation of the temperature profile Q(x). 
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